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Abstract 

As part of an effort to make use of single-mode squeezing for detection of sub-shot-noise optical 
signals, we study the balanced-heterodyne scheme, for which the corresponding spectral density of 
the photocurrent fluctuations produced at the output of the detector is calculated as the Fourier 
transform of their autocorrelation function. We show that, for maximal signal-to- noise ratio en- 
I hancement by use of squeezed states of light, an optical signal to be measured in this scheme must 

be carried in the squeezed quadrature of the carrier field. Most importantly, our analysis indi- 
cates that "the additional heterodyne noise" can be eliminated in balanced-heterodyne detection 
for single-mode squeezing under some experimental conditions and, from this, we discuss how this 
scheme may be exploited in gravitational-wave searching. To demonstrate its practical feasibility, 
we propose and study a phase-locking technique. 
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I. INTRODUCTION 



Of the essence is to sense extremely-weak optical signals in several ongoing precision 



experiments on fundamental physics 



2| and for practical applications [3|, in which the 
shot noise of light resulting from zero-point vacuum fluctuations is becoming a major limiting 
factor at the quantum level in further improvement of sensitivity. To suppress the shot noise, 
one can resort to the exotic properties of the squeezed states or the entangled states of light 
^, |5|. As a matter of fact, large degree of squeezing of more than lOdB has been achieved 
experimentally 6|. Depending on the nature of the optical beam to be measured, detection of 
optical squeezing may be carried out with different schemes. Direct detection is the simplest 
one and suited for bright amplitude-squeezed light, whereas homodyne detection scheme 
for quadrature- squeezed light. Another scheme using a cavity to rotate the quadrature 



orientation of optical sidebands in phase space relative to that of the carrier 



8| can be used 



to measure the quadrature squeezing of bright light. The practical problem with these three 
schemes is the intensity noise of light, which is classical and can be substantially suppressed 
in balanced- homodyne detection |9| . The great similarity of these schemes is their relying on 
the detection of the beats between a strong optical mode, the local oscillator or the carrier, 
and noise sidebands that are correlated with each other and symmetrically located around 
the strong mode in frequency spectrum. 

Although all the aforesaid detection schemes are playing crucial roles in experimental 
study for quantum information processing, their applications to precision measurements for 
detecting sub-shot-noise optical signals are restrained. This is because the beat note be- 
tween a slowly-varying sub-shot-noise optical signal and the local oscillator is a signal close 
to DC and will be inevitably contaminated by the low-frequency dark-current noise of a 
detector. Actually, to circumvent the low-frequency electronic noise, optical squeezing is 
usually measured at frequencies off the carrier's optical frequency in homodyne detection 
6|. Apparently, this strategy does not apply to the case where optical squeezing must be 
measured together with an optical signal that is being under investigation. To shift the 
frequency of the beat note beyond the low-frequency noise regime of the detector, one may 
consider the traditional heterodyne scheme with single local oscillator, which unfortunately 
does not work for detection of squeezed light as proven, for example, by CoUett et al. in 



10]. Heterodyne detection schemes for gravitational- wave observation 
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4l4| have been 



investigated, showing that an additional quantum-noise contribution may exist due to vac- 
uum fluctuations in frequency bands that are twice the heterodyne frequency away from the 
carrier. In the present paper, we investigate a balanced-heterodyne scheme for detecting 
sub-shot-noise optical signals, in which the "additional heterodyne noise" disappears in cer- 
tain experimental circumstances, as one will see in the following. One should note that the 
scenario here is different from that of the theoretical work for two-mode squeezing detection 



m 



15|, which does not deal with the detection of single- mode squeezing. 



The balanced-heterodyne scheme manifests itself by utilizing two local oscillators with 
equal strength that are symmetrically located in frequency spectrum around the optical 
signal to be measured, as depicted in Fig. [T]with detailed configurations to be described in 
the next section. Section III is devoted to calculating the spectral density of the photocur- 
rent fiuctuations produced by the detector. Unsurprisingly, our calculations show that the 
detected noise level of the optical signal is phase sensitive, i.e., the noise level depends on 
some relative phase of light, a feature also shared by homodyne schemes. For some given 
relative phase, observation of a sub-shot-noise signal is possible if the optical signal is carried 
in the squeezed quadrature of the field. On the other hand, if the signal is carried in the 
anti-squeezed quadrature, it must be manipulated to be carried in the squeezed quadrature, 
using for example a sing le-ended cavity QQ, beto. ,ete.od.„e detection . pe*.n,ea. 
Moreover, we discuss how the additional heterodyne noise can be ruled out in balanced het- 
erodyne detection. What succeeds in Section IV is of how to lock the relative phase of light 
for balanced heterodyning. The locking technique is quantitatively analyzed to further show 
the practical feasibility of the studied scheme. Finally, we discuss the potential application 
of the balanced-heterodyne scheme for practical gravitational-wave observation, followed by 
a short summary as conclusions to this work by the end with acknowledgements. 



II. HETERODYNE DETECTION SCHEME 

According to the uncertainty principle, when one quadrature of an optical field is 
squeezed, its canonical conjugate should be ant i- squeezed. In the conventional homodyne 
detection scheme, one may choose to measure either the squeezed quadrature or its conju- 
gate by controlling the phase of the local oscillator relative to some reference defined by the 
measured optical mode. Changing this relative phase actually orients the direction of the os- 
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cillator along that of the detected quadrature in phase space (Fig. [2^). Only when the local 
oscillator is parallel to the squeezed quadrature, could one observe the most noise-reduction 
below the shot-noise level in the homodyne detection. In contrast, the heterodyne scheme 
with single local oscillator has proven unsuitable for squeezing detection 10|], because the 
relative phase and, hence, the relative orientation in phase space, of the local oscillator with 
respect to the squeezed quadrature continuously varies with time and therefore is out of 
experimental control. As a remedy to this, the balanced-heterodyne scheme takes advantage 
of two local oscillators with equal strength, the superposition of which oscillates but remains 
pointing in a fixed direction in phase space, selecting the corresponding quadrature to be 
measured (Fig. Wp)- To see this, we may write the quantity of light intensity at the output 
port of the beamsplitter in the following form (see Fig. [T] for reference) 

= a^a + 2S cos (fit + 50)X(0) + 4(f ^ cos^ (fit + 50). (1) 

Here we adopted units in which the light intensity is expressed in photons per second. 
E(+)(t) = ae-^"^"*, 4^\^) = tfe-*('^«+^)*+'^<^i, and 4^\t) = c?e-'('^«-^)*+^<^2 represent the 
optical field to be detected and the two oscillators, respectively. And we assumed that the 
two oscillators are strongly-excited coherent modes with the same amplitude {£' a. positive 
real number) and can be approximated by classical fields. X{(j)) = X coscj) + P sin cj), where 
X = ae"*^ -I- a'l'e*^ and P = ae"*'^''"'"'^/^^ -|- d'^e^'^^'^^^'^'^ are two quadrature amplitudes of the 
optical field E^-\t). = (0i + 02)/2 - /3, and (3 is an arbitrary phase that is used to define 



X and P, thereby the major axis of t 
in the direction of the X or P axis 
the two oscillators. 



le distribution ellipse of X and P in phase space points 
6(p = (02 — 0i)/2 is the phase difference between 



A. Photocurrent fluctuations 

The first term on the right-hand side of Eq.([l]) is the DC signal of the detected field, which 
is weak compared to other terms and, hence, negligible in strong-oscillator approximation. 
The last term contributes a classical signal that does not fluctuate when the oscillators 
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are in highly coherent states 19|. Accordingly, the middle term dominates in the photon 
fluctuations: 

Ai{t) ^ 2^ cos {nt + 50)AX(0), (2) 

which shows phase sensitivity on (p and differs from the homodyne case by only a modula- 
tion function cos {Qt + 6(f)). In other words, through controlling the relative phase 0, one 
can choose to measure the noise of any quadrature of the detected field in the balanced- 
heterodyne scheme, which is analogous to the conventional homodyne scheme. Since the 
modulation function cos {Qt + 6(j)) in Eq. leaves unchanged the relative orientation of 
the measured quadrature of the detected field and the superposition of the two oscillators 
(Fig. Wp), its only effects are to split the spectral density of the photocurrent fluctuations 
into two parts and to shift each part along the frequency axis by +Q or —Q, as will be shown 
in the following. 



B. Photoelectrical signal 

At the output of the detector, the photoelectrical signal, with photon fluctuations ne- 
glected, is proportional to the light intensity averaged over the initial states of the optical 
fields. Using Eq. ([1]), we may calculate the averaged light intensity as 

< i{t) > = <d^a> +2^ cos {nt + 50) < X(0) > +4^2 cos^ {nt + 6(j)). (3) 

Apparently, the heterodyne signal of interest at frequency Q has an amplitude 2<f < X(0) > 
that exhibits a dependence on 0, in other words, a photoelectrical signal appearing in < 
X(0) > for a given may not show up in the conjugate < X(0 + tt/2) >. 

Consequently, a naive picture that one may imagine is as follows: To take advantage of 
optical squeezing to achieve maximal enhancement of signal-to-noise ratio in precision mea- 
surements, one must find a 0, for which the photoelectrical signal at heterodyne frequency 
Q is maximized and the photon fluctuations are suppressed the most at the same time. To 
put it another way, for sub-shot-noise detection of an optical signal in balanced-heterodyne 
scheme, the quadrature X(0) in which the optical signal is carried must be a squeezed one. 
True or not, we expect to find the answer in the subsequent section after calculating the 
spectral density of the photocurrent fluctuations produced at the output of the detector. 
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III. SPECTRAL DENSITY OF PHOTOCURRENT FLUCTUATIONS 



We will follow the approach of Ou, Hong and Mandel (l8| to calculate the spectral density 
of the photocurrent fluctuations. Nonetheless, stationary photocurrents can be assumed only 
at the time scale of Q~^, beyond which the measured spectral density should be treated as 



an average over measurement time T (see Ref. 10|, l20| for similar treatments). Usually 
T ~ Qr^-^ » (Qr the detection bandwidth. To resolve the heterodyne beat note 
at Q, one must set Qr << ^^). In consequence, the spectral density of the photocurrent 
fluctuations AJ{t) = J{t)— < J{t) > reads 



1 + 

7f dt dre'"^^ < AJ{t)AJ{t + r) >, (4) 

TJ g J _^ 



where < AJ(t)AJ(t + r) > is the auto-correlation function of the photocurrent fluctuations. 
Supposing that every photoelectron emitted at time t' gives rise to a definite photoelectrical 
current pulse j{t — t') for t > t', then the total photocurrent is 

i 

wherein the sum is taken over the various random-emission times ti. When t < ti, j{t—ti) = 
because no photoelectrical pulses exist yet prior to the emission of electrons. With Eq. ([5]), 



the auto-correlation function of the photocurrents can be computed as [18|, l21j 

< J(t) J(t + r) > = ^ < j(t - tl)j{t + T-tra)> 



Lm 



^<3it- U)3it + T-ti)> + ^<3{t- ti)j{t + T-tra)> 



I 

+ Y^3it- ti)jit + T- tm)P2{ti, tJAtiAtm, (6) 

where Pi(t)At and P2{t, t')AtAt' are respectively the probability of photodetection registered 
at time t within time interval At and joint probability that two photodetections are enrolled 
at t within At and at t' within At', respectively. The first term is attributable to the 
shot noise of the photoelectrical currents, while the second term depends on the fluctuation 
nature of the optical fields at the photoreceiver. These probabilities are related to the light 



intensity as |22l. |23| 



Pi{t)M = ri< i{t) > At 
P2{t, t')AtAt' = rf <3r : I{t)f{t') :> AtAt', (7) 

in which rj is a parameter characterizing the response of the detector to the incident hght. 
The symbol 3^ :: stands for time- and normal-ordering of the field operators. Converting 
the summation into integrand in Eq. ([6]) and using Eq. ([7]), one arrives at 

POO 

< AJ{t)AJ{t + T)> = r] dt' <i{t- t') > j{t')j{t' + r) 

n POO 

+rf j j dt'dt"\{t -t',r + t'- t")j{t')j{t"), (8) 

where introduced is the correlation function of light-intensity fluctuations 

X{t,L) =< £r : Ai{t)Ai{t + l) :> . (9) 

One should recall that < I{t) > in heterodyne detection is not time- independent, in contrast 
to the homodyne case where < i{t) > can be assumed to be stationary. The calculation pro- 
cedure of the intensity correlation function A(t, t) is trivial, and the treatment is resemblant 



to that based on the approach of Ou, Hong and Mandel pj|. The key idea is to establish a 
connection between the spectral density of the photocurrent fluctuations and the normally 
ordered, time-ordered correlation functions of the quadrature fluctuations of the optical field 
detected. 

For the heterodyne configuration as in Fig. [H the correlation function of intensity fluc- 
tuations is found to be different from that for the traditional homodyne configuration (see 
the appendix for detailed calculations): 

A(t, l) = < ^ : Ai{t)Ai{t + i) :> 

_j_ r(2,0)(^^-) ^^inL+i{<f>i+<f>2) _|_ g-iQt+i(</>i+<^2) _j_ ^-in{2t+L)+2i()>i _|_ gin(2t+t)+2i</.2j _|_ g j _|_ 

(10) 

wherein T^'^'^\l) =< AE^-\t)AE^+\t + l) > e''^°' and r(2'0)(i) =< AE'^-\t)AE'^-\t + l) > 
g-ja;o(2t+i)^ Eq. (ITOl) obviously contains terms that vary with a temporal period of (2fi)~^. 
Because the observed spectral density of the photocurrent fluctuations are averaged over a 



time period T oc » Q~^, these time-dependent terms will be ruled out of Eq. (jl]) by 
the temporal integrand and not show up in the spectral density. For this reason, we will 
keep only the time-independent terms in Eq. fllOl) and replace the correlation function X{t, t) 
with a t-independent function 

A'(r) ^ cos (fir){r(i'i)(r) + r(2.o)(^)e^(<^i+^2) _^ ^^^ j^ (i^) 

in the succeeding calculations. Here the strong-oscillator approximation is utilized such that 
the terms in S''^ in Eq. ffTOl) dominate over all others, which have henceforth been discarded. 
Using Eq. ([3]), similar treatment can be applied the < i{t) > term on the right-hand side of 
Eq. dHD, resulting in ^ dt < i{t) >= 2^2_ Plugging Eqs. (E]) and into Eq. (HD with 
these manipulations leads to 



= / dre 



oo p poo 



2cff'r] dt'j{t')j{t' + T)+r]' dt'dt"j{t')j{t")X'{T + t'-t") 

Jo J Jo 

(12) 

In comparison to the spectral density of the photocurrent fluctuations in conventional ho- 
modyne scheme 18|, |2l|, Eq. f lT^ differs by a global factor of two and a cos {^t) function 
(see Eq. fill I) ) in the second term. A factor of two means that the spectral density in our 
heterodyne scheme is 3dB higher, which is obviously due to the usage of two local oscillators. 
As for the function cos {^t) , as one will see, it plays a crucial role in the heterodyne scheme 
in splitting squeezing spectrum into two parts and shifting their centers off the carrier fre- 
quency, leading to the "traditional heterodyne noise" for certain experimental parameters. 

The spectral density given by Eq. ( JT2l) can be easily related to the time-ordered, normally 
ordered second-order correlation functions of the quadrature fluctuations of the detected field 
through [21 1 

Rer(i'i)(r) = [rn(r) + r22(r)]/4 
Imr(i'i)(r) = [ri2(r)-r2i(r)]/4 
Rer(2-0)(r)e2^^ = [rn(r)-r22(r)]/4 

Imr(2.o)(r)e2^^ = - [T^2{T) + r2i(r)] /4, (13) 

where r„„(r) =< 3^ : AE^(t)A£'„(t r) :>,(m,n = 1,2), and Ei{t), E2{t) are the 
quadrature operators defined as 

Ei{t) = i(+)(t)e^(^"*-^) + E(-)(t)e-^(^°*-'') 

E2{t) = ^(+)(t)e^(-<'*-/3--/2) _^ ^(-)(^)g-i(a;ot-/3-V2)_ (^4) 



[5 is the arbitrary phase associated with the field quadrature -'(t), and AE" = E— < E >. 
After substituting Eqs. f lT3|) into Eq. (fTTj) . one achieves 

A'(r) (I'^^os (r]r){rn(r)(l + cos 20) + r22(r)(l - cos 20) + [T^^ir) + r2i(r)] sin 20}, 

(15) 

= (01 + 02) /2 — /3 again. Let K{(jj) be the Fourier transform of the photoelectrical current 
pulse j(t), 

noo 

K{uj) = / rfrj(r)e*"^ (16) 

which may be interpreted as the frequency response of the detector, and let $m.n(i^) the 
Fourier transform of the correlation functions T^niT) of the field quadrature fiuctuations: 

/ + 00 
dTTmn{r)e'^^{m,n = l,2). (17) 
-oo 

With the help of Eqs. f lT^ . f ITB]) . and ffT7|) . one can readily rewrite the spectral density of 
the photocurrent fiuctuations Eq. fll2p as 

x{uj)^V'^'\K{uj)\'{2 + 

+ (r7/2) [$n(u; + ^]) + $ii(u; - r^)] (1 + cos 20) 
+ (r7/2) [$22(w + fi) + $22(w - r^)] (1 - cos 20) 

+ (r//2) [$12(0; + ^]) + $21(0; + r]) + $12(0; + $21(0; - r])] sin 20}. (18) 
Particularly, in the special case that (p = kn [k any integer), 

Xioo) ^ 27]<^^\K{uj)\^{l + {r]/2) [^n{oo + fi) + $ii(a; - Q)]}. (19) 
It follows that, if $11 < 0, the measured spectral density x{^) is to fall below the vacuum 



level for = kn, which is similar to the case of the conventional homodyne detection 18l. |21|. 
When (j) = kn ± n /2 {k any integer), 

Xicu) ^ 2ric^^\K{co)\^{l + {ri/2) [$22(w + n) + $22(0; - n)]}. (20) 

Again, akin to the homodyne scheme, the spectral density will be lower than the vacuum 
level for = A;7r±7r/2, if $22 < 0. Whether $11 or $22 is to be observed is conditioned on the 
specific values of 0, confirming the intuitive conclusion drawn upon Eq. ([2]) together with 
Eq. ([3]): To achieve sub-shot-noise detection of an optical signal with the best signal-to-noise 
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ratio enhancement in the balanced-heterodyne scheme, the optical signal must appear in the 
carrier's quadrature < X(0) > that is squeezed. 

Let one consider an optical signal carried in the amplitude of a field. The heterodyne 
signal, proportional to < X{(j)) > according to Eq. ([3]), at frequency Q is maximized 
when (j) = kn {k integer) because it is carried in the amplitude of the field. In this case, 
noise reduction below the shot noise level is available in the photocurrent only if $11 < 0, 
as is shown by Eq. ( fT9l) . meaning quadrature-amplitude squeezing. On the contrary, if 
it is the quadrature phase of the incident field that is squeezed, i.e., $22 < 0, one cannot 
observe amplitude-maximized signal with noise reduction, since maximized heterodyne signal 
demands (p = kir and, on the other hand, to observe $22 < requires = kn ± tt/2. 
Consequently, in the case of $22 < 0, before the balanced-heterodyne scheme is applied 
to yield a sub-shot-noise signal, one needs to rotate the axis of the squeezed^ quadrature 
by an angle of 90° in phase space, for instance, with a single-ended cavity [16], such that 
a quadrature-phase squeezing $22 < is converted into a quadrature-amplitude squeezing 
$11 < 0. 

On the other hand, Eqs. ( |T9l) and (l20l) both entail two equally-split parts of = 1,2), 
each shifting away from the original center by ±Q, as a direct consequence of the existence 
of cos{Qt) in Eq. f|TT]) . which is nothing but our enthusiastic argument previously made on 
Eq. ([2]). This will cause an issue in practice when applying balanced-heterodyne scheme in 
precision measurements: The spectrum of squeezing $11 (cj) or $22(1^) can put a limit to the 
degree of noise reduction in balanced-heterodyne detection. For a conceptual illustration, let 
one consider a squeezed-vacuum field generated from an optical parametric down-converter 
in a cavity. Suppose that this squeezed field is to be heterodyned with two local oscillators 
in the configuration as depicted in Fig. [H The spectrum of squeezing at the output of the 
cavity is ^] 

= ~^(7/2+7)2 + c.2' (21) 

where 7 = (1 — M)c/l is the cavity mirror refiectivity, c the speed of light in vacuum, 
and / the cavity length) is the cavity damping rate and e is a measure of the effective 
pump intensity of the down-converter. Obviously, 7/2e is a key parameter that determines 
the spectrum of squeezing. To achieve sub-shot-noise detection with signal-to-noise ratio 
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improved the most, one must control the relative phase = kir, so that the spectral density 



x{u)^2r]c^'\K{u)\' 



e7 €7 



(22) 



(7/2 + e)2 + + fi)2 (^/2 + ey + {io-n)\ 

Numerical simulations based on this equation (Fig. E]) show that, to fully take advantage 
of the degree of squeezing of the incident signal, the heterodyne frequency must satisfy 
Q << '-f. This is the condition under which "the additional heterodyne noise" can be 
completely removed in the studied scheme, at least, in principle. This may find interesting 
applications in precision measurements, such as in the observation of gravitational waves. 
In the limit of Q » 7, one can obtain a noise reduction by 3dB at most, because of 
the apparent fact that $11 (u; + Q) makes no contribution to noise squeezing at frequencies 
where $11(0; — Q) does its best, the physics behind which is the vacuum fluctuations twice 
ihe heterodyne frequency away from the carrier, an important subject of previous works 
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IV. PHASE LOCKING TECHNIQUE 

To enforce the balanced-heterodyne scheme for sub-shot-noise detection, the optical signal 
to be measured should appear in the squeezed quadrature of the carrier field, as discussed 
earlier. Meanwhile, the relative phase (p = (0i + 02)/2 — P must be well controlled at 
some values for maximal < X(0) > according to Eq. ([3]), requiring an electrical locking 



technique for this scheme. A well known technique of quantum noise locking [2^ proposed for 
homodyne scheme apparently is not applicable to the balanced-heterodyne scheme entailing 
two local oscillators. We propose to phase-modulate the two oscillators simultaneously with 
the same modulator and utilize the interference of the detected optical signal with two of the 
four sidebands of the oscillators, see Fig. HI The beat notes of the detected signal and those 
modulation-created sidebands are to be demodulated, then low-pass-filtered by a loop-filter 
to obtain an error signal for the locking loop. This phase-locking scheme is actually a version 
of the coherent-modulation-locking technique, whose stability is much better than the noise 



locking scheme for homodyne detection 24 1. 



At the output of the photodetector, the photocurrent is [21 1 

/+00 
dt'Pi{t')j{t-t') 
-00 
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/ r 1 
dt' 7] < i{t - t') > j{t') 

^7]q< i{t) >, (23) 

in which q = dt'j(t') is the total electrical charge delivered by the current pulse resulting 
from one photoelectron. Here one assumed that the electrical current pulse is much shorter 
than the oscillation period of < i(t) >, in another word, < /(t) > varies slowly compared to 
electrical pulse j(t), which is usually the case in practice. Under this assumption, < I(t) > 
may be considered approximately constant as long as the electrical pulse lasts and does not 
contribute to the integral in Eq. ( l23l) . 

After passing the same electro-optical modulator, the oscillator fields become 

wherein < ^2 is the phase-modulation frequency and 6 the modulation depth. Then 

< Jit) >^m< \E^'\-t) + 4'\-t) + 4'\-t)] X W^^\-t) + + > 

= r]q< (a^ + ^e*™-^'^^^*) + ^e''^'"''^^^'^) x (a + ^6-^™+^'^^^*) + ^e^™+^'^^(*)) > 

+ m^JM > e'^'e-'^^-^'^' - m^MO) < «^ > e'-^^e'^^-^')* + (25) 

where an appropriate modulation depth is assumed so that e*^*^™^'* ^ Jo{d) + Ji(6)e^^'^ — 
Ji{e)e-'^'' with Jo,i(^) the zero- and the first-order Bessel functions. In the last step, 
explicitly written out are the terms oscillating at frequency Vt — Vt' and the reason is that 
a demodulation signal at this frequency is to be used for error-signal pick up. Hence, the 
error signal to be demodulated reads 

r 9 - - 1 

S{(P)=Ar]qSJi{d)-^<X{(P)> sinm - n')t + 5(p\. (26) 
\-d(j) -I 

The term ^ < X(0) > in this equation ensures that (j) can be readily locked at any values 
for which < X(0) > is maximal, i.e., the peak of the photoelectrical signal. Nevertheless, 
one must note that the relative phase 50 of the two oscillators should also be kept steady 
for the locking scheme to work well. 
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V. BALANCED HETERODYNING FOR GRAVITATIONAL- WAVE OBSERVA- 
TION 



The existence of the additional heterodyne noise makes the previously-investigated 
schemes 
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14j less competitive than homodyne schemes in practical applications for pre- 



cision measurements. However, our analysis on the balanced-heterodyne scheme based on 
the quantum theory of optical coherence shows the possibility of eliminating this additional 
noise under certain experimental circumstance. Therefore, the balanced-heterodyne scheme 
may be a promising choice for photoelectrical readout in gravitational-wave observation 
experiments. The price that one pays is that the heterodyne frequency Q is not free to 
choose. Instead, the choice of fl must guarantee that the degree of squeezing at optical 
frequencies uq ± 2Q is more or less the same as at uq (Fig. |3]). 

Since the balanced-heterodyne scheme makes use of dual local oscillators with equal 
strength, the unbalanced Schnupp sidebands generated by phase modulation in Advanced 
LIGO IJ] cannot be used as the balanced local-oscillators. However, one might pick up part 



of the light from the bright port of the interferometer to serve as oscillators for balanced 
heterodyning after some appropriate manipulations. If this is the choice, the proposed phase- 
locking technique in the preceding section may be exploited to control the fluctuations of 
the involved phases of light. 



VI. CONCLUSIONS 

Although homodyne detection scheme has been widely used to measure the quantum 
noise of light below the shot-noise level in the field of quantum information processing, it can 
hardly be exploited to detect slowly- varying sub-shot-noise optical signals due to the intrinsic 
electrical noise of the photodetector at low frequencies. To fulfill the need of noise reduction 
in precision experiments where the shot noise of light is becoming a substantial restrictive 
factor for further sensitivity improvement, we have investigated the balanced- heterodyne de- 
tection scheme, where observation of noise reduction below the shot-noise level is sensitive to 
some relative phase defined by the signal- carrying optical field and the oscillators, a trait akin 
to the classical homodyne schemes. An optical signal carried in the squeezed quadrature of 
the carrier field can be measured with the best signal-to-noise ratio enhancement, whereas 
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that carried by the quadrature conjugate to the squeezed one may be manipulated by a 
single-ended cavity before a balanced-heterodyne detection is performed. Most importantly, 
we have shown that, under certain experimental circumstances, the additional heterodyne 
noise existing in previously-studied schemes may disappear in balanced-heterodyne detec- 
tion, which may be particularly interesting in applications for gravitational-wave searching. 
A phase-locking technique has been analyzed to yonder demonstrate the feasibility of the 
scheme for practical implementation. 



Appendix 



According to the definition, the correlation function of light-intensity fluctuations reads 



x{t, i) = < sr : AI{t)Ai{t + l) :> 

= < ^ : i{t)i{t + l):> - < i{t) >< i{t + l) >, 



(A.l) 



where J(t) = [E^-\t) + <st\t) + <St\t)] x [E^+\t) + ^/^^(t) + 4'^\t)] in conformity with 



Ml 



.(+)/ 



Eq. ([T]). One may expand the two terms on the right-hand side of Eq. flA.ip separately as 
follows: 

< ^ : i{t)i{t + i) :> 

= < [^(-) (t) + (t) + 4-^ (t)] X [e(-) (t + + 4'^^ (t + + 4'^ + 0] 



(+ 



X [E 



+ [4^ 

+ 
+ 

+ wt 
+ wt 
+ wt 
+ [4^ 
+ wt 
+ [4^ 



{t + L) + 4^\t + L)+ 4^\t + L)] X [EW(^) + 4'^\t) + 4'^\t)] 



> 



it) + 4 
(t) 



62 



it) 
it) 
(t) 



(- 



(- 



it) 
(t) 



4 

i$2 



it) + ^2 



(- 



(+ 



m4^ 
it)][4^ 

it)Wt 



it)][4^ 
{t)][4^ 



(t) + 4'^\t)wr\t + + '^f\t + L)][sr\t + + + 0] 

(t + + 4^"^ + 0] [^/^^ (t + + 4^^^ it + 0] < (t) > 

(t + + + 0] [^/^^ (t + + it + 0] < ^^^^ it) > 
(t) + 4^\t)][4^\t + l) + 4^\t + 1)] < E^-\t + i)> 
(t) + 4^\t)][4-\t + l) + 4-\t + 0] < E^^\t + i)> 

i)][4^\t + l) + 4^\t + 0] < E^-\t)E^^\t) > 

(t) + 4^\t)] < E^'\t + i)E^^\t + > 

(t + + 4'\t + 0] < E^'\t)E^+\t + l)> 
(t + + 4^\t + l)] < E^~\t + L)E^+\t) > 
{t + i) + 4^\t + 0] < E'^'\t)E^-\t + i)> 
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+ [4~\t + + 4'\t + 0] < E^'\t)E^+\t + L)E^+\t) > 

+ [4^\t + + 4^\t + 6)] < E^-\t)E^-\t + > 

+ < E^-\t)E^-\t + + i)E^^\t) >, (A.2) 



< i{t) >< i{t + L) > 

= < [M-\t) + 4-\t) + 4-\t)] X + + > 

X < \E'^-\t + + 4'\t + + 4'\t + l)] X + 6) + 4^\t + + 4^\t + 0] > 



= 

+ [4^ 
+ [4~ 
+ 

+ wt 
+ wt 
+ wt 
+ wi"- 
+ wt 

+ [4^ 

+ wt 
+ wt 
+ wi"- 
+ [4~ 

+ Wi^"- 

+ < E 



(t) + 4- 



(t) + 4 
it) 



4 



it) + 



(- 



(t + + 



(t) 


+ 4- 


(t 


(t) 


+ 4^ 


(t 


(t) 


+ 4- 


{t 


it) 


+ 4^ 


{t 


(t) 


+ 4- 


it 


(t) 


+ 4- 


(t 


(t) 


+ 4^ 





(t + + 4' 



(t + + 



(+ 



[4^ 
[4- 

wt 

[4^ 

[4' 

{t + 

[4^ 
[4- 

[4^ 
[4^ 

[4- 

< E 

< E 
it + 
{t + 



{t)E'^+\t) ><E 



t) + 4^\t)]wt\t + + 4'\t + ^)][4^\t + + 4^\t + 0] 
t + i) + (t + l)] [4^^ {t + L) + 4+) {t + i)] < (t) > 



t + L) + 



(t + 0] [4''^ {t+L) + 4^^ {t + 0] < {t) > 



t > 



t) + 4^\t)][4^\t + + 4^\t + 1)] < M-\t + t)> 
t) + 4^\t)][4'\t + l) + 4-\t + i)] < E^^\t 
)][4^\t + i) + 4^\t + 0] < E^-\t)E^^\t) > 

t) + 4^\t)] < E^'\t + L)E^^\t + l)> 

(-) 



t 

t + L 



L) + 6"^ '{t + L)] < E^'\t) >< E^+\t 



4 it + 0] < E'^'Kt + l) >< M+\t) > 

t + L) + 4^\t + 0] < E'^~\t) >< E^~\t + l)> 

t + L) + 4'\t + 0] < E^^\t + >< E^^\t) > 

>< E^-\t + L)E'^+\t + i) > 
-\t) >< E^-\t + L)E^+\t + i) > 
)] < E^-\t)E^+\t) >< iW(t + l)> 
)] < M-\t)E^+\t) >< E^-\t + l)> 
\t + i)E^^\t + i) > . 



(A.3) 



Since the seven leading terms in the expansion of < 3^ : I{t)I{t + l) :> are identical to those 
in < I{t) >< i{t + t) >, they are canceled out when Eqs. (IA.2p and (IA.3P are plugged into 
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Eq. (1A.1|) . Therefore, 



A(t, 6) = [4^\t) + 4^\t)][4~\t + + 4~\t + i)] < AE^-\t)AE^+\t + l)> 
+ [4'\t) + 4'\t)][4^\t + + 4^\t + l)] < AE^-\t + L)AE^+\t) > 
+ [4^\t) + 4'^\t)][4^\t + + 4^\t + i)] < AE^-\t)AE^-\t + l)> 
+ [4'\t) + 4'\t)][4'\t + + 4~\t + l)] < AM+\t + OA^W(t) > 
+ O(^), (A.4) 



where 0((f) represents the terms in ^ and those without S' in Eqs. flA.2p and flA.Sp . with 



(g' being the field amphtude of the two oscillators. Substituting 4^\'t) = ^e-^('^«+^)*+^'^i, 
and 4^\''^) = (fe-^('^«-")*+^'^2 j^to Eq. ([Oil , one can arrive at Eq. after some simple 



mathematical manipulations. 
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FIG. 1. (color online) Schematic for heterodyne detection of a sub-shot-noise optical signal. The 
optical field E^^^{t) = ae"*^"* to be detected is combined with two local oscillators, £'i^\t) = 
^^-i(u;o+n)t+i<i,i S'^'^\t) = S'e-'-('^o-n)t+i<i>2^ ^t a beamsplitter with near 100% transmissivity. 

= (S'^^\t) +S'2~^\t) stands for the superposition of the two oscillators. Even if the oscillators 
are attenuated when reflected by the beamsplitter, we still assume that they are much stronger with 
equal amplitude than the detected field. In the frequency domain, the oscillators at wq ± are 
symmetrically located around the detected field at ujq. 01^2 are the global phases of the oscillators 
and f3 is determined by the squeezed quadrature of the detected optical mode. BS: beamsplitter. 
D: photodetector. SA: spectral analyzer. 
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FIG. 2. (color online) Illustration of the relative orientation of the squeezed quadrature of the 
detected field and the local oscillator(s) that can be experimentally controlled by some relative 
phase of light, (a) Conventional homodyne scheme, where controlling the relative phase, (f)io — /3, 
of the oscillator S'e^'^'" with respect to the reference determined by the squeezed quadrature of 
the field a selects which quadrature to be measured. The axis of the squeezed quadrature is 
parallel to the direction of the oscillator in phase space whenever <j)io — 13 = kn (k any integer), 
(b) Heterodyne scheme with dual local oscillators (S'^~^\t) and (S'2^\t), in which the superposition 
f^(+){t) = (^i^\t) + 4^\t) of the two oscillators remains in a line with a fixed azimuthal angle 
relative to the squeezed quadrature of the detected field E{t). Through controlling the relative 
phase ^ = {(t>i+(l)2)/2-p defined by (+)(t) and E{t), one may select which quadrature of E{t) to be 
measured. The direction of the superposed field is parallel to the axis of the measured quadrature, 
even though each of the oscillators rotates at frequency clockwise or counter-clockwise with 
respect to the detected field. See Fig. [T]for the definition of (j)i^2 and (3. 
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co/y 

FIG. 3. (color online) Numerical simulations for noise reduction in balanced-heterodyne detection 
(using Eq. (j22p ) and homodyne detection of sub-shot-noise optical signals, assuming 7/2e = 1 for 
simplicity. Normalized x{^) refers to x('^)/(2??fo'^[-f^(w)|^), see Eq. (j22p for reference, (a) Het- 
erodyning with r2/7 = 0.05. Almost perfect squeezing is present in the photocurrent fluctuations, 
indicating the absence of the additional heterodyne noise present in other heterodyne schemes [3~ 
\l\ due to vacuum fluctuations in frequency bands that are twice the heterodyne frequency Q, away 
from the carrier, (b) Heterodyning with Vt/^ = 0.5. The degree of squeezing shows a tendency 
of degradation due to increased quantum noise at ojq ± 20. (c) Heterodyning with Vi/^ = 5. The 
greatest degree of squeezing is only 3dB (50% noise reduction), showing vacuum fluctuations at 
ujQ it 20. (d) Homodyne detection with the same squeezed- vacuum mode detected as in (a)-(c). 
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FIG. 4. (color online) Phase locking scheme for balanced-heterodyne detection of sub-shot-noise 
optical signals. The two oscillators are phase-modulated by the same modulator at a RF frequency 
0,' < Q with a modulation depth 9. The error signal is fed into a piezo-electrical transducer 
to change the global phase of the oscillators, after the photoelectrical signal is demodulated at 
frequency 17 — Q'. Obviously, the two RF signals at and il' must be phase-locked to each 
other. The control part for the relative phase of the two oscillators is omitted in the figure, albeit 
also important as explained in the text. EOM: electro-optical modulator. PZT: piezo-electrical 
transducer. BS: beamsplitter. D: photodetector. DC: electrical directional coupler. SA: spectral 
analyzer. LP: loop filter. A: amplifiers. and il': signal generators. PLL: electrical phase-locking 
loop. 
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